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Abstract 

We consider a relativistic plasma containing charged chiral fermions in an exter- 
nal magnetic field, e.g a chirally symmetric quark-gluon plasma created in relativis- 
tic heavy ion collisions. We show that triangle anomalies imply the existence of a 
new type of collective gapless excitation in this system that stems from the coupling 
between the density waves of the electric and chiral charges; we call it "the Chiral 
Magnetic Wave" (CMW). The CMW exists even in a neutral plasma, i.e. in the 
absence of the axial and vector chemical potentials. We demonstrate the existence 
of CMW and study its properties using three different approaches: i) relativistic 
magnetohydrodynamics; ii) dimensional reduction to (1 + 1) Sine-Gordon model, 
appropriate in a strong magnetic field; and iii) holographic QCD (Sakai-Sugimoto 
model), appropriate at strong coupling. We also briefly discuss the phenomenolog- 
ical implications of the CMW for heavy ion collisions. 
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1 Introduction 



Recently, the role of triangle anomalies in the dynamics of relativistic plasmas in magnetic 
field and/or at finite angular momentum has excited considerable attention. Such plasmas 
are created for example in relativistic heavy ion collisions at RHIC and LHC where the 
initial energy density significantly exceeds the threshold for the production of decofined 
and chirally symmetric quark-gluon plasma, and the coherent electromagnetic fields of 
colliding ions create a pulse of very intense magnetic field. Of particular interest are the 
following two phenomena caused in the quark-gluon plasma by the axial anomaly: the 
Chiral Magnetic Effect (CME) and the Chiral Separation Effect (CSE). 

The CME is the phenomenon of electric charge separation along the axis of the applied 
magnetic field in the presence of fluctuating topological charge [H El El HI [5] . The CME 
in QCD coupled to electromagnetism assumes a chirality asymmetry between left- and 
right-handed quarks, parametrized by an axial chemical potential /i^. Such an asymmetry 
can arise if there is an asymmetry between the topology-changing transitions early in 
the heavy ion collision. In particular, at finite axial chemical potential fi^, an external 
magnetic field induces the vector current ji = ip^ip: 

jv = -^f^AB] (1.1) 

in our present convention the current of electric charge is ejv- Closely related phenomena 
have been discussed in the physics of primordial electroweak plasma [6] and quantum 
wires [7]. While the original derivation used the weak coupling methods, the origin of 
the effect is essentially topological and so the CME is not renormalized even at strong 
coupling, as was shown by the holographic methods [H [9l [101 IHl 1121 ES] . The evidence 
for the CME has been found in lattice QCD coupled to electromagnetism, both within 
the quenched approximation ITSl 116] and with light domain wall fermions [IT] . 

Recently, STAR [THl [19] and PHENIX [201 [21] Collaborations at Relativistic Heavy 
Ion Collider reported experimental observation of charge asymmetry fluctuations. While 
the interpretation of the observed effect is still under intense discussion, the fluctuations 
in charge asymmetry have been predicted [1] to occur in heavy ion collisions due to the 
CME. Additional tests include the correlation between the electric and baryon charge 
asymmetries [22] . There is an active ongoing discussion of the microscopic mechanisms of 
CME [231 [211 [251 [261 [271 [281 [29] and of the quantitative estimates of the expected charge 
asymmetries and of possible backgrounds - see e.g. [SQl [311 1321 ESI EH [351 EEl [371 138] . 
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The Chiral Separation Effect (CSE) refers to tlie separation of cliiral cliarge along tlie 
axis of external magnetic field at finite density of vector cliarge (e.g. at finite baryon 
number density) [391 HHl SI] ■ The resulting axial current is given by 

JA = ^^^vB, (1.2) 

where /xy is the vector chemical potential. The close connection between CME and CSE 
can be established for example by the method of dimensional reduction appropriate in 
the case of a strong magnetic field |12]: the simple relations Jy = J\, = Jy between 
the vector Jy and axial Ja currents in the dimensionally reduced (1 + 1) theory imply 
that the density of baryon charge must induce the axial current, and the density of axial 
charge must induce the current of electric charge (CME); see also Ref.^3]. Since in the 
strong coupling, short mean free path, regime the plasma represents a fluid (for a recent 
review, see |11]), a number of recent studies initiated by pH!] address the effects of triangle 
anomalies in hydrodynamics, e.g. [IHl HGJ STJ HE]. 

The central observation of the present paper is the following: the connection between 
the CME and CSE implies the existence of a new type of a collective excitation in the 
plasma. This excitation stems from the coupling between the density waves of electric and 
chiral charge. Let us illustrate this statement by a qualitative argument, to be followed by 
more rigorous derivations in sections O |3l and HI Consider a local fluctuation of electric 
charge density; according to eq. fll.2p it will induce a local fluctuation of axial current. 
This fluctuation of axial current would in turn induce a local fluctuation of the axial 
chemical potential, and thus according to eq.f ll.ip a fluctuation of electric current. The 
resulting fluctuation of electric charge density completes the cycle leading to the excitation 
that combines the density waves of electric and chiral charges; we will call it the "chiral 
magnetic wave" (CMW). 

Apart from being interesting in its own right, the existence of CMW has important 
implications for the phenomenology of heavy ion collisions. The CME relies on the fluctu- 
ation of the axial charge density and so the net effect is expected to vanish when averaged 
over many events; one thus relies on measuring the fluctuations of charge asymmetries 
[U US]. On the other hand, since the quark-gluon plasma produced in heavy ion colli- 
sions possesses non-zero value of the baryon chemical potential, the CSE can lead to a 
non-vanishing axial current even after the summation over events is performed. However 
unfortunately a direct detection of the axial current in heavy ion collisions is very chal- 
lenging. The CMW should exist even in a neutral plasma and so can induce interesting 
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observable effects in heavy ion collisions even after the sum over many events is performed; 
we will return to this topic in the Summary. 

The paper is organized as follows. In section |2] we provide a derivation of the CMW 
based on relativistic magnetohydrodynamics. In section |3] we consider the case of a strong 
magnetic field and perform a dimensional reduction; in this case the dynamics of CMW 
is described by the Sine-Gordon equation. In section H] we describe the CMW at strong 
coupling using the holographic methods within the Sakai-Sugimoto model. When the 
electromagnetism is treated dynamically, the CMW mixes with the longitudinal charge 
wave in the plasma - the plasmon. We consider the mixing of CMW with plasmons in 
section [6l Finally, in the Summary we outline the main result of the paper and discuss 
the directions for future studies. 

2 Chiral magnetic wave in magnetohydrodynamics 

Let us now proceed with the derivation sketched out in the introduction. We will see that 
there indeed exists a new gapless excitation in a deconfined QCD plasma that propagates 
along the apphed magnetic field; it arises as a dynamical consequence of the underlying 
triangle anomaly of chiral symmetry. This new excitation is a long wavelength hydrody- 
namic mode with a dispersion relation that looks like that of sound waves, 

uj = ^v^k - iDik'^ + ■■■ ; (2.3) 

however, these propagating modes carry both electric and chiral charges. Since these 
modes would not exist if it were not for the applied magnetic field or the underlying 
triangle anomaly, we will call them "the chiral magnetic waves" (CMW). They give rise 
to several important new transport properties of hot QCD plasma, and affect its ther- 
modynamics; we will further discuss this in section [61 We note that CMW exists even if 
the background plasma is neutral under either baryonic or axial symmetry, which should 
make it a generic phenomenon in relativistic plasmas. 

For simplicity, let us consider single flavor (Np = 1) massless QCD with chiral sym- 
metry U{1)l X U{1)r, or equivalently U{l)v x U{1)a where V{A) denotes vector(axial) 
respectively. The axial symmetry U{1)a suffers from both QCD anomaly with gluonic 
topological density and from the triangle anomaly of global chiral symmetry. The latter 
is in fact not harmful to the conservation of U{1)a as long as one does not elevate the 
global chiral symmetry to a gauged one, while the former indeed breaks the axial U{1)a 
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symmetry by quantum fluctuations of topological density. 

Our starting point is tlie anomalous generation of vector and axial currents along the 
applied magnetic field in the presence of axial (vector) chemical potential /ia (/iy), as 
given by eqs. (11 .ip and (ll.2p . We will now re- write these equations in a more suggestive 
matrix form as 

We are interested in small linearized fluctuations of the chiral currents Ja and jv in the 
plasma; let us assume that this plasma is neutral, without any background charge density 
on average. We may then perform a linear expansion of the chemical potentials with 
respect to small charge densities (jy, j^), 




J \ W W I \ 3a J V ^" W V aAV ^AA I V A 



+0 ((/)') 



(2.5) 



Remembering that 



~ Q„-nQ„-n ■ (2-7) 



37 

/^^ = 7wJ , ^ = y.A (2.6) 

where T is the Helmholtz free energy, the a's appearing above are nothing but the sus- 
ceptibility matrices of vector/axial charge densities, 

Considering the parity V transformation V — )■ — V and A — )■ A, one concludes that parity 
invariance of QCD implies that ay a = (^av = in the neutral plasma, [ly = \ia = 0. 
Moreover, a simple large iVc counting shows that 

ayy ^aAA-O (^^^ , (2.8) 

while their difference in a deconfined and chirally symmetric phase is subleading 

aw -OiAA--0 ■ (2.9) 

we will confirm this within the holographic large Nc Sakai-Sugimoto model in section |H 
Independently of this, the vanishing of the difference ayy — aAA can be taken as a signal 
of chiral symmetry restoration. Therefore, we expect it to be a good approximation to 
let ayy = aAA = a in the chirally symmetric phase; this leads us to 



jy\^ KeBa / l\(3'y 
JA 2vr2 \1 J{ j\ 



(2.10) 
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It is natural to diagonalize the equation above by going to the chiral basis 

fL^lu^-jl) , A^lu^ + fA) ■ (2.11) 

In terms of chiral currents, our previous assumptions and the definition of a's are easily 
translated to 

1 /(9/ii\ 1 f d'^T \ 1 fdi2R\ 1 / d'^T 



The fl2.10p then leads to two decoupled relations 

jL,R = Ti-^^jjlR , (2.13) 

where one should keep in mind the definite sign in front of the right-hand side depending 
on the chirality of the currents. 

One can view the above expression as the leading constitutive equation for the currents 
in the long wavelength derivative expansion of hydrodynamics. Indeed, our starting point 
(]l.l|1.2p is strictly valid only when the variation of chemical potentials is sufficiently 
slow; for a finite frequency/momentum these expression gets modified resulting in fre- 
quency/momentum dependent chiral magnetic conductivity [501 El EH |32] . The equation 
fl2.13p is the first leading term in the derivative expansion, while the next leading-order cor- 
rection to the chiral magnetic conductivity will be or u'^ ~ k"^ in frequency/momentum 
space. However, there is an important first-order derivative term in any constitutive equa- 
tion of conserved current: a diffusion term — DVj°, with a diffusion constant D. In our 
case, we will be interested only in the waves propagating along the magnetic field direction 
which we call longitudinal; thus on general grounds, the constitutive relation including 
the next leading-order diffusion term reads as 

3l,r = t\-^^\31r-Dl ' jIr + --- , (2.14) 

with a longitudinal diffusion constant Dl. Although we discuss only longitudinal dynam- 
ics in this paper, it would also be interesting to study the transverse dynamics with the 
transverse diffusion constant Dt- 

A similar constitutive equation was written previously by Son and Surowka [4]J. There 
is however one point that will appear important for us: while Ref.|j41J considers a weak 
magnetic field B and treats it in the linear approximation, we are claiming that f l2.13p 
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and (I2.14P are valid for arbitrary strength of eB non-perturbatively. This is equivalent 
to the validity of our starting point f ll.ip . fll.2p for arbitrarily large eB, which is not at 
all trivial and is a consequence of the absence of corrections to the axial anomaly. It is 
also important to note that although fll.ip .( !L2|) look linear in eB, this linearity is only 
apparent. Given a fixed density jy^, the chemical potentials fiv,A in general may well 
depend on the dynamics of underlying microscopic theory, such as coupling constants, 
temperature, as well as magnetic field eB non-linearly, so that the currents in (ll.ip . fll.2p 
can in fact be very non-linear in these parameters, see e.g. [511 [53]. The statement 
of (11. ip .( [1^21) is that these dependencies can be absorbed into the chemical potentials 
liv,A- Therefore, one expects that a and Dl are in general non-linear functions of eB, 
temperature T, etc. In particular, they would also depend on the coupling constant, 
so that it is meaningful to study them in the strong coupling regime via holographic 
QCD as we do in our section [H Let us mention Ref. [5^ that presents a diagrammatic 
proof of (ll.ip . (ll.2p perturbatively in the coupling constant, and Ref.[T^ that proved this 
relation in the Sakai-Sugimoto model using the two-derivative approximation. We will 
present a strong coupling proof with full DBI action of holographic QCD in section [H with 
arbitrary strength of eB, which presumably includes non-linear effects of derivatives as 
well. Therefore we expect that the equations (ll.ip . (ll.2p and hence (12.140 hold universally. 

Our next step is to combine (I2.14p with the conservation law d^j'^ = 0. We take 
B = Bx^ and consider only longitudinal gradient di, which results in 

{d,T^^^d,-D,di)jl^ = ^ . (2.15) 

This describes a directional wave, or chiral wave, of charge densities whose direction of 
motion is correlated with its chirality. The velocity is given by 

_ N.eBa _ IUB_ ( d^xiA _ 
'''' 27r2 47r2 [dfj An^ [df^J " ^ ^ ^ 

As we discussed in the previous paragraph, one expects that and are non-trivial 
functions of eB, T, and the coupling constant, so they are interesting dynamical quantities 
to compute in any model. In frequency/momentum space, the above equation takes the 
form 

uj = zpv^k - iDLk"^ + ■ ■ ■ , (2.17) 

as a hydrodynamic dispersion relation. Our main observation is the new first term in 
the dispersion relation which makes the mode propagating instead of simply diffusing. It 
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exists only if 1) triangle anomaly exists and 2) there is a background magnetic field. We 
stress that this new chiral mode of electric and chiral charge transport is present even if 
the plasma is neutral on average. 



3 Chiral magnetic wave in strong magnetic field: 
dimensional reduction to (1 + 1) Sine-Gordon prob- 
lem 

The expression f l2.16p for the velocity of the Chiral Magnetic Wave (CMW) shows that 
~ eB. What happens when eB becomes large? We will now show that in the limit 
eB —7- oo the velocity stays finite and reaches the velocity of light. To understand 
this, let us first examine the spectrum of charged fermions in magnetic field; for massless 
fermions, the energies of Landau levels are given by 



Er,= \l2eB{n + ]^-s-,)+pl- (3.18) 

we assume that magnetic field B is directed along the z axis. Note that the lowest Landau 
level (LLL) with zero energy (at vanishing momentum p^) is not spin-degenerate, and all 
excited Landau levels are - therefore the net chirality Sz ■ Pz is carried only by the LLLs. 

When magnetic field is large compared to the temperature T, eB ^ T^, the fermions 
stay "frozen" in the LLL and the transverse and longitudinal (along z) dynamics are inde- 
pendent. In this case one can perform a dimensional reduction to the (1 + 1) dimensional 
theory where the only allowed direction of motion is along the magnetic field. The trans- 
verse density of states is given by eB/(2'K), the density of LLs in the transverse plane. 
The longitudinal phase space density for a Fermi-momentum is simply pf/{2tt). There- 
fore for the density of Nc left-handed fermions j° with pf = f^i (for massless particles the 
Fermi-momentum and the chemical potential /z are equal) is 

fL = N."^^, (3.19) 



and the derivative is given by 



djl eB 
Substituting this expression into fl2.16p . we obtain 



' (3.20) 



= 1; (3.21) 
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therefore in the hmit of strong magnetic field eB ^ the CMW indeed propagates with 
the velocity of light. 

We can deduce a more detailed information about the dynamics of the CMW in strong 
magnetic field by making use of bosonization procedure. As is well known, bosonization 
approach is very powerful in the studies of (1 + 1) dimensional systems [55| [56] . This is 
easy to understand as in one spatial dimension the produced fermion and an anti-fermion 
never separate and propagate together as a composite bosonic excitation - even if they 
do not interact at all! 

The conservation of vector current dfj^jy = can be ensured by introducing a boson 
field and choosing 

= -^e^'^d,^; (3.22) 



this way the vector current is always conserved independently of the equations of motion. 
The corresponding choice for the axial current is 

fA = e'^J^ = -^d^^. (3.23) 
\/7r 



For zero quark mass and in the absence of background electric field, the axial current 
should be conserved: 

d.j'^ = n^ = 0. (3.24) 

Therefore the conservation of axial current leads to the wave equation for the bosonic 
excitation ip. The variation of in space and time causes variations of both charge and 
chiral densities; therefore it is natural to identify the wave defined by (13.241) with the 
CMW. Since (for massless quarks) there is no mass term in (13.241) . this CMW propagates 
with the velocity of light. 

It is convenient to split the field cp into left- and right-moving components: 

(p{x,t) = (Pl{x - vt) + (pji{x + vt); (3.25) 
in terms of these fields the original fermion fields are 

^^(x - vt) = ^= e^^v^v^^C^-^-t). ^^(3; + = J_ ^^Vi^^L{x+vt) _ ^3 25) 

v27r v27r 



Upon quantization of the boson fields, the fermions (I3.26P represent coherent states of 
Vr^^l- Viceversa, these boson fields describe the collective electric and chiral charge 
density fluctuations of the underlying fermion quark fields. 
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(a) 



(b) 



Figure 1: Anomalous triangle diagrams of chiral magnetic waves for (a) abelian and (b) 
non-abelian flavor symmetries in the presence of external magnetic field B. The quark 
lines are dressed propagators including the background magnetic field. 

The quark mass term in the original Lagrangian can be re-written with the help of 
fl3.26p in the following form: 



and leads to the interaction term for bosons; the bosonic Lagrangian thus describes the 
Sine-Gordon theory. This correspondence between the fermionic and bosonic description 
has been discovered by Coleman [55] and Mandelstam [56] . The kinks of the Sine-Gordon 
theory describe the fermions, and the fluctuations of the boson field - collective fermion- 
anti-fermion excitations. Coleman's theorem [57] forbids spontaneous breaking of a con- 
tinuous symmetry in (1 + 1) dimensions. The underlying reason for this no-go theorem 
is the presence of strong infra-red fluctuations that dominate two-point functions and 
destroy the long-range order in two dimensions. This means that the massless (in the 
limit of massless quarks) boson is not a Goldstone bosoij*]. This is consistent with our 
interpretation of the field as of a collective density wave of electric and chiral charges. 
At finite density of baryon or chiral charge, this density wave propagates on top of a "chi- 
ral spiral" [59] - the winding configuration of the background field if causing the chiral 
magnetic effect in (1 + 1) dimensional description [^ H3] . 

*There exists however a way around Coleman's theorem uncovered by Witten |58) : for A^-component 
field at iV — cx), when the number of degrees of freedom diverges at each point, the two-point functions 
can exhibit long-range order and the Goldstone phenomenon can still be realized. 




(3.27) 
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One can also generalize the above to non-abelian version of chiral magnetic waves. 
As FigurelU illustrates, essentially the same type of triangle diagrams for mixed U{l)v x 
SU{Np)v X SU{Nf)a would result in 



N,eBa f 1 \ f fy^ 



(3.28) 



for non-abelian SU{Np) components of currents jyj^ = iy^jst°- with tipit^'t^) = 
One then gets to the same conclusion on the emergence of the chiral/directional CMWs 
for each non-abelian component of SU{Nf)- Upon the 1 + 1 dimensional reduction with 
strong magnetic field, these non-abelian chiral magnetic waves should be described by 
non-abelian bosonization [60] of SU{Nf) symmetry. More precisely, the theory of Np 
Dirac fermions in fundamental representation of SU {Nc) in bosonized description can be 
represented as 

C = N^C{SU{Nf)) + C{U{1)) + NfC{SU{N^)) , (3.29) 

where C{G) represents the Wess-Zumino-Witten model of group G with level 1 [60]. After 
integrating over QCD SU{Nc) dynamics, one is left with the first two pieces as a low 
energy effective theory [611 |62l [631 El]- The U{1) part is what we have discussed above, 
while the non-abelian SU{Np) part describes the non-abelian chiral magnetic wave as a 
propagating group field g{x) G SU{Nf). 



4 Chiral magnetic wave in holographic QCD 

In this section, we intend to study chiral magnetic wave at strong coupling in the frame- 
work of holographic QCD, particularly using the model by Sakai and Sugimoto [65]. The 
Sakai-Sugimoto model is a right place to look for the phenomenon because it includes 
the relevant triangle anomalies of QCD chiral symmetry in terms of 5D Chern-Simons 
terms and the existence of chiral magnetic wave is robust as long as the right anomalies 
exist. What is non-trivial will be the details of dispersion relations such as wave velocity 
and diffusion constant, which do depend on the strong dynamics and the stength of the 
applied magnetic field. 

The Sakai-Sugimoto model in deconfined phase consists of two separated D8 and D8 
branes that touch the black hole horizon of the background geometry as shown in Figure 
[21 This geometric separation of two probe branes correctly indicates that interactions be- 
tween left-handed quarks and right-handed quarks are sub-leading in large expansion. 
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X_4 



Black-hole horizon 



Figure 2: A schematic picture of D8—D8 branes in deconfined phase of the Sakai-Sugimoto 
model. 



when there is no chiral condensate in the deconfined phase. Because of this decoupling, 
one can study dynamics of each D8 probe brane independently at least in leading order 
approximation. The background geometry is a warped product of a 5D black hole and 
X where DS/D8 branes are separated along S^. They wrap all other dimensions 
except to make up 9 dimensional world-volume on which the action is given by a 
Dirac-Born-Infeld action plus a Chern-Simons term; 

SDs/m=-l^s J d'C e~f Vdet {g* + 27r/gF) t ^ gf J F.^^'aAaFAF , (4.30) 

with Hp — {27r)~Pl7^'^^^ ■ After integrating over S"^, one arrives at an effective 5 dimen- 
sional world- volume action of DS/D8 branes embedded in the 5 dimensional black hole 
space-time of the metric 

(^)' (|-/(C/)rfi' + E(rf.xOj +2dC/di , /([/) = !- (^^y , (4.31) 

where U is the holographic direction whose UV(IR) boundary is sitting at C/ — > oo{U — 
Ut) , and the metric is written in Eddington-Finkelstein coordinate. Note that Eddington- 
Finkelstein coordinate is suitable to describe future event horizon in a smooth manner, 
and the in-going boundary conditions of perturbed modes at the future horizon reduce to 
simple regularity in this coordinate. The temperature T of the dual 4 dimensional QCD 
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plasma is related to the parameters in the above by 



The resulting 5 dimensional world-volume action of the probe brane is 
5^g/-gg = -Ci?t j d^xdU U-^ yjdet {gl^ + 27r/2F) 

1 

T d'^dUe^^^^^^^AuF^pFQu , C = ,(4.33) 



where the sign of the last Chern-Simons term depends on the chirality of the probe D8/D8 
branes. The e-symbol above is numerical one. It is straightforward to write down the 
equation of motion of the world-volume gauge field Am from the above action, which 
reads as 



cmnild, 



•n 



f/Vdet {g* + 2nP,F)(^{g* + 27illF)-') 



(4.34) 



where [MA^] = MN — NM is anti-symmetrization. In deriving the above as well as for 
later convenience, it is useful to be reminded of the following expansion 

Vdet (1 + M) = 1 + ^tr {6A) + I [tr {6A)]' - hi ({6Af) + O ({6A)') . (4.35) 

2 8 4 

For our purpose of studying chiral magnetic wave, we first need to construct a back- 
ground solution of having a constant magnetic field of arbitrary strength pointing, say, 
direction. It is in fact easy to show that the trivially constant F12 = B solves both Bianchi 
identity and the equation of motion f l4.34p . so the holographic background configuration 
of a constant B will be our starting point of looking for a signal of chiral magnetic wave. 
In making contact between our D8/D8 branes and the QCD chiral symmetry, we work in 
the case of Np = 1 for simplicity, neglecting f/(l)yi-anomaly induced by gluons as a lead- 
ing large approximation. This is appropriate as long as we focus on triagle anomalies 
which are encoded in the 5D Chern-Simons terms in holographic QCD. Therefore, chiral 
symmetry is U{1)l x U{1)r, and D8(L'8)-brane dynamics captures U{l)L{U{l)f() chiral 
dynamics of QCD holographically. The electromagnetism lies in the diagonal combination 
of U{1)l and f/(l)/j, while axial symmetry is the other orthogonal combination, so one 
has a dictionary 

cAem = li^L + Ar) , Aa = ^ {-Al + Ar) , (4.36) 
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where AiiAfj) is the external potential that couples to the chiral current Jl{Jr), and e 
is the electromagnetic coupling constant. Therefore, having a constant electromagnetic 
B means having a constant F12 on each D8 and D8 brane with 

F,f = Fj^=eB , (4.37) 

and we will assume this in the following. 

Our next task is to expand linearly around this background to study hydrodynam- 
ics of longitudinal charge/current fluctuations. As we discuss in previous sections, these 
longitudinal charge/current fluctuations would have had leading diffusive dispersion re- 
lation u ~ —iDk"^ in the absence of triangle anomalies. Chiral magnetic wave is an 
anomaly-induced modification of this into a leading propagating dispersion relation 

u: ~ ^v^k - iDLk'^ + ■■■ , (4.38) 

where the velocity as well as the longitudinal diffusion constant are expected to 
depend on strong coupling dynamics and the magnitude of the magnetic field B. The 
above dispersion relation due to triangle anomalies is our main objective we will be heading 
to in this section. Especially, as a function of B will be a new physical quantity that we 
compute in the framework of holographic QCD. Because should vanish in the absence 
of anomalies and the magnetic field, it will be an odd function of the product of Chern- 
Simons coefficient and the magnetic field B. We will confirm this expectation later. We 
also emphasize that the sign of the leading term in (14.381) is fixed by the sign of Chern- 
Simons term, that is, by the chirality one is looking at, so left-handed charge/current 
fluctuations would have a definite propagating direction with respect to the magnetic field 
direction, while the right-handed ones would propagate in the opposite way. Therefore, 
one can achieve chirality separation in this way as a dynamical consequence of triangle 
anomalies of QCD. Recall that no charge in the background, either axial or baryonic, needs 
to be present for the effect to take place, and only magnetic field applied to a deconfined 
plasma with an underlying triangle anomaly is sufficient to have the phenomenon. 

Studying linearized fluctuations from the constant F12 = eB background requires 
expanding the 5D action (I4.33P quadratically in terms of relavant fluctuation fields, and 
it is tedious but straightforward to do it using the formula (I4.35p . As we are interested 
in longitudinal charge/current fluctuations, it is sufficient to consider {6Ftu,SFsu,SFt3) 
fluctuations of "helicity" only. Because of residual 5*0(2) rotation symmetry of the 
constant F12 background, other non-zero helicity modes simply decouple from the above 
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modes at linearized level of equations of motion. After a sizable amount of computation, 
one arrives at 



o(2) _ 

D8/D8 



d^xdU - [A{U) {6Ftuf - B{U) {6F,uf + 2C{U) {6F,u)] 



j d^xdU [SAuSFts - SAsSFtu + SASFs, 
where the three functions that appear in the coefficients are given by 

A{U) = C {27rllY U + {2nlleBy 
B{U) = C {2711^^)^ Uf{U) + (2<e5)' 



(4.39) 



C{U) = C(27r/2)V 



(4.40) 



The second line in f l4.39p is from the Chern-Simons term that represents triangle anomaly 
of QCD holographically. It is easy to keep track of that by a combination N^eB which 
is the product of anomaly coefficient and the magnetic field. From the above quadratic 
expansion of the action, one easily writes down the linearized equations of motion as 



eB 

du {A{U)5Ftu) + C{U) {d,6F,u) T ^^'^^3C/ 

eB 

du {C{U)6Ft3) T -^SFtu 



du iB{U)6Fsu) + C{U) idt6Fsu 
A{U) (dtSFtu) - B{U) {d,6F,u) + C{U) {d,6Frs) ± 



47r2 
NccB 



6F 



t3 









(4.41) 



To proceed to hydrodynamic analysis from the above, one invokes low frequency/momentum 
expansion in solving fl4.4ip order by order with right boundary conditions. As mentioned 
before, simple regualarity at the future horizon f/ = f/^ in Eddington-Finkelstein coordi- 
nate is equivalent to implementing in-coming boundary conditions, and one should also 
impose normalizability on the modes at the UV boundary [/ — > oo. Near f/ — )■ oo, fl4.4ip 

3 

gives two asymptotic behaviors, either 0{1) or (9(f/~2), and the normalizability picks 
up the latter only. Assuming the frequency/ momentum factor e~*'^*+*'^^^ or equivalently 
replacing dt = —iu and = ik, and working in a gauge Au = 0, the equations of motion 
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become 

N eB 

du [A{U) {du6At))] + zkC{U) {du6A,) t (duSA,) = 0, 

du [B{U) (duSA^)] - luCiU) (St/Ma) - du [C(U) {iuj5A^ + ik5At)] T {duSAt) = 0, 

eB 

-iojA{U) {du6At) - ikB{U) {du6A^) + ikC{U) {iuj5A^ + ik6At) ± {iu5A^ + ik6At) = 0, 



(4.42) 



which should be solved in perturbative expansion of {u,k). What one expects is that 
imposing boundary conditions restricts the solution space such that u in hydrodynamic 
expansion is determined once k is given, and the relation u = u{k) is called the dispersion 
relation. 

There are several methods in literature to solve similar kinds of equations in hydro- 
dynamic expansion [66], but we will follow our own method which seems most convenient 
to us. We first assume that uj{k) is an analytic power series in k, which is expected based 
on hydrodynamics, 

co{k) = ank"" = v^k - iDLk"^ + ■ • • . (4.43) 

n>l 

Inserting this to (14. 42 p . then one can take k as the only expansion parameter in solving 
(I4.42P systematically, along which a„ should also be determined order by order. Because 
(I4.42P is linear in {5 At, 5A-i) one can always rescale them so that they start their k- 
expansion as 



(M„ 5A,) = ^(aS"\ 4"))r = (Af ), ) + {Ai'\ A^^)k + ... , (4.44) 



^iW _ ^ /i(o) /,{o)^ , ^4(1) 

n>0 

where {A^f\ A^^^) cannot vanish simultaneously by definition. It is straightforward to 
insert (I4.43P and (I4.44p into (I4.42p . and solve order by order in k. 

At C(/c°), one gets the equations 



du 

du 



MU){duAr^)]T^{duAf^)=0 , (4.45) 



N^eB 

m){duAf^)]T^{duAf)=0 , (4.46) 
-^v,AiU) [duA?^) - ^B{U) [du^^) ± ^ (^v,Af + .^f )) = . (4.47) 
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Integrating the first two equations fl4.45p and fl4.46p gives us 



AiU){duAf)T'-^AP = C, , (4.48) 
BiU)lduAf^)T^A?^ = , (4.49) 

with two integration constants, while the last equation fl4.47p simply becomes 

-tv^Ci-iC2 = ^ ^^x = -^ . (4-50) 

so that will be determined once Ci,2 (or more precisely their ratio) are fixed by imposing 
relevant boundary conditions. Considering fl4.48p at f/ — )■ oo, one first fixes 

C, = lim A{U) (duAfA , (4.51) 

and solving A^^^ gives 

Inserting this into (14.49 p . one gets a second order differential equation for A^'^^ 



t 1 

2 



BiU)3„ (AiU) (a,Ml»>)) - (^) Al»> ^ ±C, (^) . (4.53) 

Because B[Ut) = at the horizon, the regularity boundary condition imples that 



C2 = T[^]Ar{UT) , (4.54) 



and one can write the solution for Ai as 



(0) 

t 

(0) _ 4 , 4(0)/ 



Ar = A + Ar{UT) , (4.55) 
where A satisfies 

B{U)du{A{U)(duA))-[^^' A = Q , (4.56) 

with the boundary condition A{Ut) = at the horizon. This uniquely determines A up 
to rescaling. Note that Af\UT) is free up to this point, and the final boundary condition 
we need to impose is to demand vanishing Af''' at U — )> oo, and from fl4.55p this fixes 

A^Vt) as 

aI°\Ut) = - lim A{U) , (4.57) 
16 



so that C2 is finally 

'N,.eB 



= . (4.58) 

Observe that Ci is also given by A as 

Ci = lim AiU) (duA) . (4.59) 

l/-5>00 \ / 

Therefore, the complete solution at 0{k!^) with the right boundary conditions can be 
written solely in terms of A as above, and it is unique up to overall rescaling. Especially 
Vy^ is well-posed and given by 

As expected, is proportional to the anomaly coefficient, and its sign depends on the 
chirality and the chiral magnetic wave is uni-directional. 

Discussions in the previous sections independently argue that should be given by 

where (/i, j°) are chemical potential and charge density for either U{1)l or U{1)r. Indeed 
we can confirm this expectation in our final formula fl4.60p which is fully non-linear in 
ei?, so that one can consider this as a strong coupling proof of the relation. First note 
that to compute j one only needs a linear perturbation of /i or j° to the system 
of our background B field, and the relevant equations of motion for them are precisely 
given by our previous one fl4.4ip with additional assumption of space-time homogeneity 
dt = ds = 0. Then the equations simplify exactly to the previous f l4.48p and f l4.49p 
with suitable boundary conditions. In the second equation considering the horizon point 
U = Uti we have B{Ut) = and we demand that At vanishes at the horizon, so that 
C2 = 0. Combining the two equations after removing ^3, one easily arrives at that At 
satisfies the same equation f l4.56p that A satisfies, and moreover they share the same 
boundary condition at the horizon U = Ut, so that they are in fact the same object 
At = A. On the other hand, gauge/gravity dictionary tells us that up to linear order, 

/i = lim AtiU) , f = lim A{U) (duAt) , (4.62) 



so that one has 
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1.5 



2.0 



eB/(GeV)^ 



Figure 3: Numerical result for in the Sakai-Sugimoto model with T = 150 MeV(dotted), 
T = 200 MeV(plain), and T = 250 MeV(dashed). 



which proves the relation fl4.6ip at strong coupling fully non-linear ly in eB. As the validity 
of (I4.6ip is equivalent to fll.ip . fll.2p . this constitutes a holographic proof of f ll.ip .f lL2|) . 

It is interesting to see how depends on the magnitude of the magnetic field eB in 
a non-linear way because {A{U),B{U)) contain eB as in f l4.40p . although the necessary 
analysis inevitably involves numerical study. To perform numeric analysis, we have to 
specify parameters of the model. First of all, one can always put 27tII = 1 for simplicity 
because this factor will eventually cancel out in any well-defined field theory observables. 
One can easily check that this is the case for as well. By fitting to the observed p-meson 
mass and the pion decay constant, Sakai-Sugimoto fixed the parameters as 

N, = 3 , (?2^,iV,~17 , MKK--OMGeV , (4.64) 

where ^kk are related to Qs by 

gYM = '^nlsMKxgs ■ (4.65) 

Recall that what matters for us is simply the parameters and C, and in terms of the 
above parameters, one has 

C ~ 0.0211 , ~ 1.44 . (4.66) 

For the temperature, we take T = (150, 200, 250) MeV as an illustrative purpose. Note 
that this model has deconfinement phase transition at Tc = ^^^^ ~ 150 MeV [67] . We 
plot our numeric result of f ^ as a function of eB in Figure [3l 
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' ' ' ' ^ eB/(GeV)2 

1.0 1.5 2.0 

Figure 4: Numerical result for Dl in the Sakai-Sugimoto model with T = 150 
MeV(dotted), T = 200 MeV(plain), and T = 250 MeV(dashed). 

Limited analytic results for are available for two extreme regions, either eB 
or eB — 7- oo. For this purpose as well as an easier numerical analysis, it is convenient to 
consider the combination 

V{U) = "^^f ^ , , (4.67) 



in terms of which the equation f l4.56p becomes a simple first order differential equation 

with a boundary condition V{Ut) = 0. From (14.601) . is then simply given by 

^x = T(^)vioo) , (4.69) 



. 47r2 ^ 

which seems technically much easier. With this formulation, it is not difficult to derive 
the following results; 

Weak field limit : eB — > 
Strong field limit : eB — )■ oo 

-)■ =f1 (speed of light) . (4.71) 

Note that the strong field limit gives us the same result that one expects in the weak 
coupling Landau level picture that is discussed in section El It might come as a surprise 
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because the dynamical degrees of freedom in our holographic model are still mesonic 
degrees of freedom represented by world-volume gauge fields on D8-branes. They are 
color singlet, hadronic degrees of freedom, while they reproduce something related to 
underlying quark degrees of freedom. This looks quite similar to an idea of qaurk-hadron 
duality, and it would be interesting to explore this further. 

One can go on to the next order expansion in k to compute the longitudinal diffusion 
constant Dl. As the computation is conceptually straightforward, we only present our 
numerical results in our Figure HI For eB = 0, it reduces to the known diffusion constant 
in the model 

while in general it is a complicated function of eB and the temperature T. As one can 
observe in Figure HJ Dl generally gets decreased as eB increases, which seems physically 
sensible because a larger magnetic field would align quasi-particles more efficiently, so that 
their microscopic longitudinal velocity diffusion would be smaller. What is interesting is 
that for some range of eB in a low temperature, say T = 150 MeV, one seems to have a 
negative value of Z^^. Looking at the dispersion relation 

uj = ^v^k - iDik'^ H , (4.73) 

this signals an instability for sufficiently small k or long wave-length fluctuations, which 
is precisely similar to Gregory-Lafiamme instability in a gravity system [68]. In the case 
of gravity system, Gubser-Mitra conjecture [69] links this dynamic instability with a 
thermodynamic instability, and it would be interesting to study whether our instability is 
also related to some kind of thermodynamic instability. We leave this for future studies. 



5 Dynamical Electromagnetism : Mixing Chiral Mag- 
netic Wave with Plasmons 

The previous sections treat electromagnetism as providing only a non-dynamical external 
magnetic field, and neglecting dynamical nature of electromagnetic field. This is valid 
in the limit e — > while keeping eB finite. However, to describe the real world more 
precisely, it would be desirable to go beyond this approximation including dynamical 
electromagnetism. In general, dynamical electromagnetic field in the plasma couples to 
longitudinal fluctuations of vector charge density inducing the plasma waves. The chiral 
magnetic wave (more precisely its projection onto the vector f/(l)y-part) also involves 
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longitudinal charge density fluctuations; it is thus natural to expect them to mix with 
each other, resulting in interesting modifications of their dispersion relations. This will 
be the main topic of this section. We also point out that for non-Abelian SU{Np) chiral 
magnetic waves there is no mixing with plasmons as far as one does not introduce gauging 
of SU{Nf), so that the previous discussion stays intact for non-Abelian CMWsH 

Having dynamical electromagnetism, one first needs to include the Maxwell equation 
as a new dynamical equation of motion in addition to the conservation laws of currents, 

- SM'^) = ej(; . (5.74) 

What is important for us is how to determine the current jy. Typical discussions of plasma 
waves proceed by assuming a linear retarded response of the plasma to the electromagnetic 
field A", 

j^ = eU^''{u,k)A, , k^U^'' = , (5.75) 

where we work in the frequency- momentum space. Combined with fl5.74p this would 
result in a self-contained equation for from which one can extract the plasmons. Let us 
review this procedure briefiy as we are going to extend it by including the chiral magnetic 
waves and diffusion terms. Assuming a definite frequency-momentum e-«'^*+«*:^^ ^ t/^q focus 
on longitudinal polarizations {Ao,Ai), the only ones where we expect the emergence of 
plasmons. Because one has a gauge freedom of shifting by k^, one can use this to 
remove Ai and work with Aq only, which simplifies the analysis significantly. From f l5.74p 
and f l5.75p . one then obtains 

eA' = e'U'>\u,k)A' , kuA' = e'll'\u,k)A' . (5.76) 

These two equations are in fact equivalent as can be seen using the Ward identity k^U^^^ = 
uU^^—kU^^ = 0, and one has a non-trivial solution of A^ only if {u, k) satisfies the plasmon 
dispersion relation, 

e = e^U^\uj,k) . (5.77) 

One expects plasma waves in long wave-length regime A; — ?■ with finite frequency, so one 
expands H^^ in powers of ^ as 

n»"^.s^;(i+c^„^;)+... , (5.78) 



^To be more precise, EM charge is a sum of U{1) and I3 of SU{2)p for A''^ = 2, so the situation can 
be more comphcated in general. We leave a full discussion as a simple extension of ours to the readers. 
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upon which the dispersion relation becomes 



= eWo + clffe + 0{k'') ■ (5.79) 

this looks like a massive excitation of mass squared = c^cJq with an effective speed of 
light Ceff. The parameter uq is typically of the order of temperature T, and the plasma 
frequency is about Up ~ eT which can be small compared to T in weak-coupling limit; the 
plasmon thus can be an important ingredient in describing the hydrodynamic properties 
of the system. 

We will have to modify this picture in two aspects for our case of a constant back- 
ground magnetic field by including chiral magnetic waves and diffusion terms. Because 
the magnetic field B is constant while the plasma waves involve only longitudinal elec- 
tric field fluctuations, we can decouple them, and treat from now on as an additional 
fluctuation of the gauge field on top of the background magnetic field. Looking at the 
chiral magnetic effect fl2.10p . one first notices that the anoma/ows/?/ induced spatial cur- 
rents are not directly related to retarded response to the gauge field fluctuations A^; they 
are induced from the total charge densities jy^^ without any regard to how these charge 
densities appear. We will discuss these charge densities below. Also, we have to modify 
(12.1 op for finite frequency-momentum by introducing chiral magnetic conductivity a{uj, k) 
[50t E] because it generally depends on {u, k), 

Jv\_^^f a{u,k) \ffv\_ N^eBa ( (^) ^ f 



(5.80) 



where ctq = |^ is the zero frequency- momentum limit we have used before. In general, 
at finite {u, k), we expect a ^ a. Therefore the total spatial currents jv,A will represent 
the sum of the retarded response to given by (15.750 and the anomalously induced 
contribution (I5.80p . 

A = en'^A'>^!^('-^)f{.-^kD,f, , (5.81) 

where we keep our focus on longitudinal components only, and we also included the 
diffusion terms proportional to D^. Note that we also include the induced axial current 
from response to A^ through H^v because this term indeed exists in the presence of 
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background magnetic field B (think of triagle diagram of external B and which couples 
to axial current). It is important to keep in mind that the charge densities jyj^ appearing 
on the right hand side are total charge densities that may come from both response to 
A^^ as well as additional fluctuations due to chiral magnetic effects. Once we write down 
f l5.8ip . we don't and can't specify charge densities because they are free up to dynamical 
equations of Maxwell equation and current conservation laws. To be more precise, we 
have three dynamic equations; u = 0-component of Maxwell equation and two current 
conservation laws for jy^ (the u = 1 Maxwell equation becomes equivalent to z/ = 
once jy-conservation is imposed due to gauge invariance). They are homogeneous linear 
equations in terms of three variables [Aq, jy, j^), so that non-zero solutions exist if and 
only if the 3x3 coefficient matrix has zero determinant. This constraint on {u, k) will 
give us the dispersion relation. 

Although it is not necessary, it is convenient to decompose jy as 

fy = eU^'A^ + 6fy , (5.82) 

to visualize additional fluctuation 6jy to the retarded response explicitly. Let us then 
write down the three independent dynamical equations mentioned above. The = 0- 
component Maxwell equation is 

fc^A" = e + e5j° , (5.83) 

while the vector current conservation, = 0, looks as 

- zcoSfy + zk^^^ (-) j° + k'D, (en°M° + 6fy) = . (5.84) 

Finally, for axial current conservation, there is an important modification to its conser- 
vation due to triangle anomaly we are considering (note that we are still neglecting QCD 
anomaly from gluons). Recall that axial current becomes anomalous in the presence of 
non-zero electromagnetic E ■ B ^ due to triangle anomaly, 

d,fA='^,^''-'F,.F^, = '^E.B . (5.85) 

Remember that our Nc quarks have charge e in this paper. We already have a background 
magnetic field B = Bx^, while dynamical longitudinal plasma fluctuations we are consid- 
ering have a longitudinal electric field fluctuation Ei = diA^ = ikA^ in Fourier space, so 
that one has locally non-vanishing E ■ B that affects axial current conservation law as in 
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fl5.85p . The resulting (modified) conservation law of gives us 

-^uf^ + ^keU%A' + ^k^^ ( ^) (en°M° + 6fy)+eD,f^ = ik^^A' , (5.86) 

where the right hand side is the anomalous contribution that we discussed. The above 
equations fl5.83p . fl5.84p . and fl5.86p are the main equations for {A^,6jy,j'^) from which 
one can obtain dispersion relations. 

As an easy application as well as an illustration, let us turn off the magnetic field for 
a moment and consider the diffusion effects only, which may be called diffusive plasmons. 
In this case, one has Il^fy = due to S = 0, and axial current decouples with the usual 
diffusion cu = —iD^k^, while {A^,6jy) system becomes 

{tuj-eDL)6j^ = A;2eDin°M° , (5.87) 

which has non-zero solutions if and only if 

{e - e'n°°(u;, k)) {tuj - ^Dl) = k^e^DLU'^^uj, k) , (5.88) 

which gives the dispersion relation. Upon expanding 11°'' as in fl5.78p . one can solve the 
above for small k as 

u^ = ul+ [clff - iUpDL) k^ + 0{k^) , (5.89) 

where Up = euo is the plasma frequency. 

Going back to our interesting case of non-zero magnetic field i? 7^ 0, it is straight- 
forward to study f l5.83p .f l57Ml) . and f l5.86p in complete generality, but we will restrict 
ourselves to the case with = for simplicity in this paper, leaving their full analysis 
including to the future. The system then becomes 

eA' = e^U'^^A^ + eSfy , 

-tu6j^ + ik^-^(-)f^ = , 
2vr^ \croJ 

-^cof^ + ^keU%A' + ^k^^^(^){eU''A' + 6fy) = ik^^A' , (5.90) 
which mixes all three fluctuations together. From the first equation, one can replace A^ 
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with Sjy, and inserting it into the other two equations, one gets 

Ztt^ \croJ 

-'-A + {- ) i ^r^^ = o . 

from which one gets the dispersion equation 



\ a a{u),k) NcBa a{u),k) AV \ ' / 



where = as before. This equation is our master equation that governs mixing 

between chiral magnetic waves and plasma waves. 

For small magnetic field B and u ^ Up ^ eT <^ T , one expects that chiral magnetic 
conductivities are approximately the zero- frequency value a ~ a ^ ctq, and moreover 
anomaly triangle diagram gives us 

U%{u,k)^^^^^ as u,k^O , (5.93) 

so that the numerator in the right hand side of (15.921) becomes simplified. One then uses 
the previous expansion (15.781) of to solve the above equation to get 



P+{< + clff)k' + 0{k') , (5.94) 



which describes effects from chiral magnetic wave to the plasma waves. Note that the effect 
exists only with finite k, and this makes sense because chiral magnetic waves disappear 
in /c — limit. 

Another interesting limit is an infinitely large 5 — )■ oo limit, where one expects effec- 
tive reduction to 1+1 dimensional theory. In fact, making electromagnetism dynamical 
corresponds to 1+1 dimensional QED with A^^ massless Dirac fermions, or the Schwinger 
model [To] . It has been known for long time that the photon in the model becomes massive 
due to 1+1 dimensional axial anomaly, 

ml = ^ , 5.95 

' 71 

where e^j^j- is an effective 1+1 dimensional QED coupling constant. As the 4-dimensional 
triangle anomaly (I5.85P correctly reduces to 1+1 dimensional axial anomaly in the pres- 
ence of background magnetic field B, one should be able to reproduce this Schwinger 
phenomenon from our master equation (15.921) in the limit -B — )■ oo. 
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One can be more quantitative to test this connection. To find Cgjj, it is useful to 
consider a transverse area of 

27r 



d'xT = — , (5.96) 
eB 

to have a single lowest Landau level system per each 4D fermion because the transverse 
density of LLL is eB/27T. Thinking of fermion kinetic term, the proper normalization 
between 4D fermion and 2D fermion is 

[eB 

^4D = \ —'ip2D , (5.97) 
V 2,71 

which will be useful shortly when we discuss about n°°. Because one is looking at only 
longitudinal dynamics of 4D U(l) gauge field, the only relevant dynamical field is Fqi, 
and the gauge field action indeed reduces to 1+1 dimensional QED action as 

^ / d'x {F,,f = ±^ J d\ {F,,f ^^fd'x (Fo.r , (5.98) 

so that eljj = and the Schwinger photon mass should be 

9 N^e^B 

< = ^ ■ (5.99) 

To reproduce this from our equation (15.921) . note that as i? — )■ oo, the expected uj"^ = 
is also infinite and one naturally expects that chiral magnetic conductivities go to zero 
in this limit as the system cannot respond to arbitrary fast perturbations; (cr, a) — )■ as 
cu — 7- oo. Therefore, the solution of (I5.92p in this limit is found simply by demanding that 
the denominator vanishes or 

e = e^Y]P\uj,k) , (5.100) 

where should be given by the effective 1+1 dimensional theory. To find it, recall from 
(15:971) that j^^ = If JaD' so that 



As we have Nc 1+1 dimensional Dirac fermions, IIj?) is Nc times that of a single Dirac 
fermion, which can be found most easily by bosonization to a single real scalar field 
such that 




(5.102) 
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where is normalized to have a standard kinetic term L = ^d^cpd^^cj). Then, Ilg?) is easily 
computed as 

^2D = — ^10 ^10 = 2 = 2 ' 5.103 

so that the equation flS.lOOp becomes 

AT P^n i-2 

which indeed gives us u'^ = + k'^ reproducing the Schwinger model result. 

Therefore the plasmon in the dimensionally reduced theory can be seen as a result of 
the interaction of the dynamical photon with the chiral magnetic waves. 



6 Summary 

We have demonstrated that the Chiral Magnetic and Chiral Separation Effects (CME 
and CSE) in relativistic plasmas subjected to magnetic field imply the existence of a 
new type of a collective excitation in the plasma - the Chiral Magnetic Wave (CMW). 
This excitation represents the density waves of electric and chiral charge coupled by the 
triangle anomaly. In strong magnetic field the CMW propagates with the velocity of light, 
— 7- 1. In weak magnetic field, the velocity decreases; the result of the holographic 
computation is shown in Figure [31 At weak coupling, this decrease of the velocity of the 
CMW can be understood as originating from the admixture of the excited Landau levels. 

The existence of CMW in the quark-gluon plasma has important implications for the 
phenomenology of heavy ion collisions. The CME relies on the fluctuation of the axial 
charge density and so the net effect is expected to vanish when averaged over many events; 
one thus relies on measuring the fluctuations of charge asymmetries [U HQ]. On the other 
hand, the CMW should exist even in a neutral plasma, and so does not require the presence 
of the axial or baryon chemical potentials. Since it represents the coupled density waves 
of electric and chiral charges propagating along the direction of the applied magnetic field 
(that in heavy ion collisions is perpendicular to the reaction plane), the CMW can induce 
dynamical, reaction plane dependent, fluctuations of electric charge. The azimuthal angle 
dependence of these fluctuations will be determined by the wavelength of the CMW 
excitation. We will return to the consideration of phenomenology related to CMW in 
heavy ion collisions in a forthcoming publication j71] . 
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